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Abstract. Embedded clusters usually contain young stellar objects in different evolutionary stages. We investigate number 
ratios of objects in these classes in the star-forming regions p Ophiuchi, Serpens, Taurus, Chamaeleon I, NGC 7129, IC 1396A, 
and IC 348. They are compared to the temporal evolution of young stars in numerical simulations of gravoturbulent fragmenta¬ 
tion in order to constrain the models and to possibly determine the evolutionary stage of the clusters. It turns out that Serpens is 
the youngest, and IC 348 the most evolved cluster, although the time when the observations are best represented varies strongly 
depending on the model. Furthermore, we find an inverse correlation of the star formation efficiency (SFE) of the models with 
the Mach number. However, the observational SFE values cannot be reproduced by the current isothermal models. This argues 
for models that take into account protostellar feedback processes and/or the effects of magnetic fields. 
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1. Introduction 

Almost all stars form in clusters. Embedded clusters contain 
various types of young stars, making them ideal laboratories 
to study star formation as they provide a large and genetically 
homogeneous sample (see Lada & Lada l2003l for a review). 

Four classes of young stellar objects (YSOs) are distin- 
guisbed according to tbe properties of tbeir spectral energy 
distributions (SEDs). Originally Lada ( fTWTl defined Class 1, 

2, and 3 objects according to tbe slope in tbe SED from 1 to 
20 pm. Note that Class 2 sources correspond to classical T Tauri 
stars and Class 3 objects to weak line T Tauri stars. Later tbe 
extremely embedded sources (Class 0 objects) were added to 
this classification, and their observational properties are de¬ 
fined e.g. in Andre et al. J2000t . These four classes are inter¬ 
preted as an evolutionary sequence from Class 0 to 3. There 
are, however, concerns that some Class 0 sources might be 
mimicked by Class 1 objects seen edge on (Men’shchikov & 
Henning 119971 . Furthermore there are works suggesting that 
Class 2 and Class 3 objects are actual of the same age (e.g. 
Walter fT986t . Class 0 sources are deeply embedded protostars, 
possessing a large sub-mm (T > 350 pm) to bolometric lumi¬ 
nosity ratio (Lsmin/Lboi > 0.005). The Class 0 stage is the main 
accretion phase and lasts only a few 10“* yr. Class 1 objects are 
relatively evolved protostars. They are surrounded by an accre¬ 
tion disc and a circumstellar envelope. Pre-main-sequence stars 
in Class 2 and 3 are characterised by a circumstellar disc (opti¬ 
cally thick in Class 2, optically thin in Class 3) and the lack of a 
dense circumstellar envelope. The progenitors of these forming 
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stars are prestellar cores (starless cores, prestellar condensa¬ 
tions). These are gravitationally bound, dense molecular cloud 
cores with typical stellar masses that may already be in a state 
of collapse, but have not formed a central protostellar object 
yet. 

Star formation in molecular clouds is controlled by the 
complex interplay between interstellar turbulence and self¬ 
gravity (Vazquez-Semadeni et al. 120001 Larson 120031 Mac 
Low & Klessen 120041 and references therein). The supersonic 
turbulence ubiquitously observed in Galactic molecular gas 
(Blitz [1993 generates strong density fluctuations with gravity 
taking over in the densest and most massive regions (e.g., Sasao 
119731 Hunter & Fleck Elmegreen ll993l Padoan rr^ 
Klessen l^OOll Padoan & Nordlund l2002> . We call this process 
gravoturbulent fragmentation. In a cloud core where gravita¬ 
tional attraction overwhelms all opposing forces from pressure 
gradients or magnetic fields, localised collapse will set in. The 
density increases until a protostellar objects forms in the centre 
and grows in mass via accretion from the infalling envelope. 
The gravoturbulent models of molecular cloud evolution dis¬ 
cussed here can describe the entire collapse of a cloud core 
and the build-up of a stellar cluster as a function of time (see, 
e.g., Klessen et al. Eunni Klessen IMJTl Heitsch et al. IMJTl 
Schmeja & Klessen 120041 Jappsen & Klessen I2004t . These 
models provide a “snapshot” of the cluster at any time, allow¬ 
ing the comparison with observed clusters, and possibly the 
determination of the cluster’s evolutionary status according to 
the models. This permits to constrain the models and possibly 
to determine the evolutionary stage of the star-forming region 
by comparison with the models. 
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Complete unbiased surveys of clusters for the content of all 
four evolutionary classes and prestellar cores are difficult to un¬ 
dertake, since they require different observational techniques. 
Hence, different investigations have to be combined. The result 
therefore might suffer from different detection limits or vary¬ 
ing spatial coverage. With this caveat in mind, we search the 
literature for information on the different YSO classes in sev¬ 
eral star-forming regions and tried to construct roughly homo¬ 
geneous samples of Class 0 to Class 3 sources and prestellar 
cores to compare them with our models. 

The purpose of this paper is twofold; On the one hand, we 
compile an observational sample of absolute numbers of YSOs 
belonging to the different classes for several star-forming re¬ 
gions from the literature (Sect. |3, on the other hand we anal¬ 
yse the evolution of the YSO classes in gravoturbulent models 
(Sect. 13 and compare it with the observational data in Sect.|3] 
Finally, in Sect.|3we summarise and conclude our hndings. 

2. Observational data 

The detection and classihcation of prestellar cores and 
Class 0/1/2/3 objects requires different observational tech¬ 
niques. Thus, we have to construct our samples from various 
sources. We consider corresponding areas on the sky, but the 
caveat remains that the combined samples may be far from 
complete and not homogeneous. Since it would hardly effect 
the relative numbers, the problem of unresolved binaries can 
be neglected. Besides, we cannot resolve close binaries in the 
SPH models either. The absolute and relative numbers of YSOs 
adopted from the observations for the subsequent analysis are 
listed in Table [2 Prestellar cores are hard to determine both 
from observations and in our models, and in particular they are 
hard to compare, since the status of the cores (Jeans-critical or 
subcritical) is often unknown. Some cores considered as “star¬ 
less” might even in fact turn out to harbour embedded sources 
(Young et al. l2004t . Therefore we do not use them for the actual 
comparison, but keep them as an additional test for consistency. 
Due to constraints from the models (see Sect. |3 we consider 
Class 2 and 3 combined. So in the lower panel of Table ^only 
Class 0, Class 1, and the combined set of Class 2 h- 3 objects are 
shown. 

2.1. p Ophiuchi 

The p Ophiuchi molecular cloud is the closest and probably 
best-studied star-forming region, offering the most complete 
sample of YSOs. Bontemps et al. lIMin hnd a total number 
of 16 Class 1 sources, 123 Class 2 sources, 38 Class 3 sources, 
and 39 Class 3 candidates. They did not detect the previously 
known two Class 0 objects lying within their area. Following 
the reasoning of Bontemps et al. and the hndings of 

Grosso et al. ( Eunnb that there might be almost as many Class 3 
as Class 2 objects, we add the candidates to the Class 3 sam¬ 
ple, yielding a total number of 77 Class 3 sources. Stanke et al. 
( I2004t performed a 1.2 mm dust continuum survey of the p Oph 
cloud and detected 118 starless clumps and a couple of pre¬ 
viously unknown protostars. To obtain a reasonably homoge¬ 
neous sample when combining the Bontemps et al. lIMlTt and 


Table 1. Absolute (upper panel) and relative (lower panel) 
numbers of YSOs of different classes from observations (ref¬ 
erences see text) 


Region 

prestellar 

0 

1 

2 3 

p Oph 

98 

2 

15 

111 77 

Serpens 

26 

5 

19 

18 ~20 

Taurus 

52 

3 

25 

108 72 

Chal 

<71 

2 

5 

175 

IC 348 


2 

2 

261 

NGC7129 


1 

20 

80 

IC 1396A 


2 

6 

47 1 

p Oph 


0.01 

0.07 

0.92 

Serpens 


0.08 

0.31 

0.61 

Taurus 


0.01 

0.12 

0.87 

Chal 


0.01 

0.03 

0.96 

IC 348 


0.01 

0.01 

0.98 

NCC7129 


0.01 

0.20 

0.79 

IC 1396A 


0.03 

0.11 

0.86 


Stanke et al. J2004> surveys we restrict ourselves to the range 
covered by both investigations, which is the 45' x 45' region 
of the main p Oph cloud (L1688) investigated by Bontemps et 
al. J2001> . This area contains two Class 0 sources (Froebrich 
120051 Stanke et al. l2004> . The global star formation efficiency 
(SFE) in L1688 is estimated at 6 -14%, although the local SFE 
in the subclusters where active star formation takes place is 
signihcantly higher with ~ 31% (Bontemps et al. I20()lt . The 
measured velocity dispersion is 2.6 km s * in p Oph A and 
2.7 km s * in p Oph B (Kamegai et al. l2003t . With the reported 
temperatures of 11 and 7.8 K this corresponds to Mach num¬ 
bers of Af ~ 13.5 and Ai ~ 15.5, respectively. 

2.2. Serpens 

The Serpens Cloud Core, a very active, nearby star-forming re¬ 
gion, contains 26 probable protostellar condensations (Testi & 
Sargent n'998> and hve Class 0 sources (Hurt & Barsonv ri9961 
Eroebrich I2005> . Kaas et al. detected 19 Class 1 and 

18 Class 2 objects in the central region (covering the held in¬ 
vestigated by Testi & Sargent ll998t . This is an unusually high 
Class 1/2 ratio compared to other regions. Kaas et al. J2004> 
cannot distinguish between Class 3 sources and held stars and 
are therefore unable to give numbers for Class 3. Preibisch 
(I2003t performed an XMM-Newton study of Serpens and de¬ 
tected 45 individual X-ray sources, most of them Class 2 or 
Class 3 objects. Considering this and the argumentation of 
Kaas et al. J2004t we can assume that there are at least as many 
Class 3 as Class 2 sources in the relevant area. The local SEE 
in sub-clumps is around 9%, the global SEE is estimated to be 
around 2 - 10% (Kaas et al. 120041 Olmi & Testi IM)^ . The 
measured velocity dispersion is 0.3 - 0.6 km s ’ (Olmi & Testi 
I2002> . corresponding to 1 < Af < 2.5 at T = 20 K. 
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2.3. Taurus 

The Taurus molecular cloud shows a low spatial density of 
YSOs and represents a somewhat less clustered mode of 
low-mass star formation. The numbers of Class 1, 2, and 3 
sources in Taurus are reported as 24, 108, and 72, respec¬ 
tively (Hartmann I2002t . In the same area there are 52 prestel- 
lar cores (Lee & Myers I1999t . and one Class 0 and three 
Class 0/1 sources (Froebrich I2005I I. For our analysis we di¬ 
vide the Class 0/1 objects into two Class 0 and one Class 1 
object. Since the prestellar cores in the sample of Lee & Myers 
were selected by optical extinction, their number may be un¬ 
derestimated compared to other regions. Estimates for the star 
formation efficiency vary between 2% (Mizuno et al. ll995> and 
25% in the dense filaments (Hartmann l2002t . The velocity dis¬ 
persion is 0.49 km s * (Onishi et al. I1996t . corresponding to 
Al ~ 2.5, adopting a mean temperature of ~ 11 K, i.e. a sound 
speed of 0.2 km s“^ 

2.4. Chamaeleon I 

The Chamaeleon I molecular cloud harbours 126 confirmed 
and 54 new YSO candidates, most of them classical or 
weak-line T Tauri stars (Class 2/3). Furthermore, four prob¬ 
able Class 1 protostars are detected by the DENIS survey 
(Cambresy et al. rT9^ . Persi et al. (EOOTt find two more 
Class 1 sources. One object is classified as Class 0, and one 
as Class 0/1 by Eroebrich (|2005ll. Counting the latter as Class 0 
gives a total of two Class 0 and five Class 1 sources. Haikala 
et al. j2005t detected 71 clumps (some of them associated with 
embedded protostars) and a mean line width of 0.62 km s * in 
the clumps, corresponding to Al ~ 3 at T « 11 K. The SEE in 
Cha I is about 13% (Mizuno et al. ll999> . 

2.5. 1C 348 

The young nearby cluster IC 348 in the Perseus molecu¬ 
lar cloud complex contains 288 identified cluster members 
(Luhman et al. unnii, including 23 brown dwarfs. The majority 
of the objects are believed to be in the T Tauri stage (Class 2/3) 
of pre-main sequence evolution (Preibisch & Zinnecker l2004> . 
The active star formation phase seems to be finished in the 
central parts of the cluster, but southwest of the cluster cen¬ 
tre a dense cloud core containing several embedded objects is 
found. There are two Class 1 objects (Preibisch & Zinnecker 
[200^ . one confirmed Class 0 source and one Class 0 or 1 
source (Proebrich l2005t . which we count as Class 0 in our anal¬ 
ysis. Subtracting the brown dwarfs, we adopt a number of 261 
Class 2/3, two Class 1 and two Class 0 objects. The average 
velocity dispersion is 1.04 km s ’ (Ridge et al. I2003t . corre¬ 
sponding to Al ~ 5.2 at T ss 11 K. 

2.6. NGC 7129 and IC 1396A 

With NGC 7129 and IC 1396A we include two additional star¬ 
forming regions in our analysis. Although, there is no informa¬ 
tion on prestellar cores in the literature, sufficient data about the 
population of young stellar objects can be found. The numbers 


of YSOs in these two clusters are based on data from the Spitzer 
Space Telescope, which are not available for the other regions. 
Since Spitzer is very sensitive to the earliest YSO classes, the 
obtained number ratios might be overestimated in favour of 
Class 0 and 1 compared to the other regions. 

In NGC 7129 Muzerolle et al. ll2()()4ll detected one Class 0 
object (classified as Class 0/1 by Eroebrich I2005t . 12 Class 1 
objects and 18 Class 2 sources in their Spitzer data. They miss 
the core cluster members and estimate a total of 20 Class 0/1 
and 80 Class 2 objects, which is a similar ratio as in Taurus or 
p Oph. The average velocity dispersion is 1.58 km s“’ (Ridge 
et al. l2003t . corresponding to Al ~ 8 at T ss 11 K. This region 
has also been studied by Megeath et al. J2004t . see also below. 

The Elephant Trunk Nebula, IC 1396A, was investigated by 
Reach et al. (I2004t using Spitzer data, revealing three Class 0/1, 
five Class 1, 47 Class 2, and one Class 3 object. These num¬ 
bers are similar to the numbers of YSOs found by Eroebrich et 
al. (I2005bt . who however cannot distinguish between Class 1 
and Class 2/3 sources. Eor our analysis we divide the Class 0/1 
objects into two Class 0 and one Class 1 object. Reach et al. 
(I2004t estimate a SEE of 4 - 15%. 

2.7. Other star-forming regions 

Megeath et al. ll2004t report Spitzer results of the four young 
stellar clusters Cepheus C, S171, S140, and NGC 7129. They 
find ratios of Class 1 to Class 2 objects between 0.37 and 0.57, 
which is significantly higher than in the other regions listed in 
Tabled except Serpens. That indicates that these clusters are 
very young, although the same caveat for Spitzer data as above 
applies. Since no data of Class 3 objects are available we do 
not include these clusters in our analysis. 

3. The models 

We perform numerical simulations of the fragmentation and 
collapse of turbulent, self-gravitating gas clouds and the re¬ 
sulting formation and evolution of protostars as described in 
Schmeja & Klessen J2004II . We use a code based on smoothed 
particle hydrodynamics (SPH; Monaghan ll992t in order to re¬ 
solve large density contrasts and to follow the evolution over 
a long timescale. The code includes periodic boundary condi¬ 
tions (Klessen ll997> and sink particles (Bate et al. I1995> that 
replace high-density cores while keeping track of mass and lin¬ 
ear and angular momentum. The periodic boundary conditions 
ensure that, independent of the box size, all formed stars remain 
in the simulation, also later-type objects that are known to be 
more widely distributed. We determine the resolution limit of 
our SPH calculations using the Bate & Burkert ( ITWTl criterion. 
This is sufficient for the highly nonlinear density fluctuations 
created by supersonic turbulence as confirmed by convergence 
studies with up to 10^ SPH particles. It should be noted, how¬ 
ever, that the Bate & Burkert (llWTI l criterion may not be suffi¬ 
cient for describing the evolution of linear perturbations close 
to equilibrium, as suggested by Klein et al. (|200l. 

Our simulations consist of two globally unstable models 
that contract from Gaussian initial conditions without turbu¬ 
lence and of 22 models where turbulence is maintained with 
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Fig. 1. The temporal evolution of the fractions of YSO classes for three models; M2k4 (left), M6k2a (middle), and M10k2 
(right), upper panel: absolute numbers, lower panel: relative numbers. Solid line: Class 0, dashed line; Class 1, long dashed 
line: Class 2+3, dotted line: prestellar cores (only shown in the upper panel). The abscissa gives the time in units of the global 
free-fall time. The zero point of the timescale cotTesponds to the time when gravity is “switched on”. Note the different intervals 
covered in the models. The vertical lines in the lower panel indicate the time r, when the model shows the best agreement with 
the observations of Serpens (dotted), NGC 7129 (dashed), Taurus and IC 1396A (dash-dotted), p Ophiuchi (long dashes), Cha I 
(dash dot dot dot), and IC 348 (second dotted line). Note, because of similar relative YSO numbers the best-fit evolutionary times 
for Taurus and IC 1396A are essentially indistinguishable. 


constant rms Mach numbers Ad, in the range 0.1 < M. < 
10. We distinguish between turbulence that carries its energy 
mostly on large scales, at wavenumbers 1 < k < 2, on interme¬ 
diate scales, i.e. 3 < k < 4, and on small scales with 7 < k < 8. 
The naming of the models, G1 and G2 for the Gaussian runs, 
and MAlkk (with rms Mach number M and wavenumber k) 
for the turbulent models, follows Schmeja & Klessen 
Details of the individual models are given in their Table 1. 

The dynamical behaviour of isothermal self-gravitating gas 
is scale free and depends only on the ratio a between internal 
energy and potential energy: a - Fint/lFpotl- This scaling factor 
can be interpreted as a dimensionless temperature. We convert 
to physical units by adopting a physical temperature of 11.3 K 
corresponding to an isothermal sound speed Cs = 0.2 km s ', 
and a mean molecular weight p - 2.36, corresponding to a typ¬ 
ical value in solar-metallicity Galactic molecular clouds. With 
an average number density n(H 2 ) = lO^cm^^, which is con¬ 
sistent with the typical density in the considered star-forming 
regions, the total mass in the two Gaussian models is 2311 Mq, 
and the size of the cube is 3.4pc. The turbulent models have 
a mass of 1275 Mq within a volume of (2.8pc)^. The global 
free-fall timescale is Tff = 10^ yr. If we instead focus on in¬ 
dividual dense cores like in p Oph with n(H 2 ) ~ lO^cm^^, 
the total masses in the Gaussian and in the turbulent mod¬ 
els are 231 Mq and 128 Mq, respectively, and the volumes 


are (0.34 pc)^ and (0.28 pc)^, respectively. The global free-fall 
timescale is Tff = 10^ yr. Note that the number of stars is not 
influenced by the adopted physical scaling. For further details 
on the scaling behaviour of the models see Klessen & Burkert 
(120001 and Klessen et al. ll2000t . 

The YSO classes are determined as follows; The beginning 
of Class 0 is considered the formation of the first hydrostatic 
core. This happens when the central object has a mass of about 
0.01 Mq ('Larson l2003t . The transition from Class 0 to Class 1 
is reached when the envelope mass is equal to the mass of the 
central protostar (Andre et al. l2000> . The determination of the 
end of the Class 1 stage is more difficult, since this is usu¬ 
ally done via spectral indices in the near-infrared part of the 
SED. Generally, after the Class 1 stage the objects are consid¬ 
ered classical T Tauri stars that become visible in the optical. 
Hence we determine the transition from Class 1 to Class 2 when 
the optical depth of the remaining envelope becomes unity 
at 2.2 pm (K-band). Using the evolutionary scheme of Smith 
(I2000t and the standard parameters as described in Froebrich 
et al. (12005 at the end of Class 0 corresponds to a mass of 
M* a! 0.43 Mend, where Mend denotes the final mass of the 
star. The end of Class 1 is reached when M, a; 0.85 Mend- Note 
that the exact value of the mass at the transition from one phase 
to the next does not influence our results significantly. Even a 
change of the opacity value by a factor of four results in a devi- 
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ation of the corresponding mass of a few per cent only. Lacking 
a feasible criterion to distinguish Class 2 from Class 3 objects, 
we consider both classes combined. The same is done for the 
observational data. 

Finding and defining prestellar cores is a more difficult 
task. Usually one considers roughly spherical symmetrical den¬ 
sity enhancements containing no visible traces of protostars 
(e.g. Motte et al. 119981 Johnstone et al. 12000b . We attempt 
to follow this procedure and define prestellar cores from the 
Jeans-unstable subset of all molecular cloud cores identified in 
our models. We use a three-dimensional clump-finding algo¬ 
rithm to determine the cloud structure, similar to the method of 
Williams et al. Further detail is given in Appendix A 

of Klessen & Burkert <20001 . This procedure matches many of 
the observed structures and kinematic properties of nearby star¬ 
less molecular cloud cores (see the discussions in Ballesteros- 
Paredes et al. l2003l and Klessen et al. 12005b . 

In order to avoid problems with the described mass criteria 
for the determination of the classes for low-mass objects and to 
be consistent with the observations, we only consider protostars 
with a final mass Mend >0.1 Mq, which roughly corresponds 
to the detection limits of the observations reported in the lit¬ 
erature. For the same reason we only consider models with a 
numerical resolution of at least 200 000 particles. Furthermore, 
in order to get reasonable numbers of protostars in the different 
classes we select only those models where more than 37 proto¬ 
stars with Mend >0.1 Mq are formed. This reduces our set of 
models to 16. Again, see Table 1 of Schmeja & Klessen <2004b 
for further details. 

In some models, a small fraction of protostars gets highly 
accelerated (e.g. by ejection from a multiple system). Due 
to the adopted periodic boundary conditions in our calcula¬ 
tions, these objects cross the computational domain many times 
while continuing to accrete. In reality, however, these proto¬ 
stars would have quickly left the high-density gas of the star¬ 
forming region and would not be able to gain more mass. We 
therefore consider accretion to stop after the object has crossed 
the computational box more than ten times. Varying this dis¬ 
tance does not significantly influence our conclusions, e.g. in¬ 
creasing or decreasing it by a factor of two changes the num¬ 
bers derived in the next Section by less than 1%. 


4. Discussion 

4.1. The evolutionary sequence 

Figure [0 shows the temporal evolution of the fractions of the 
different YSO classes for three selected models. The formation 
of the entire cluster takes place on varying timescales between 
about two and 25 global free-fall times. In the Gaussian col¬ 
lapse models and the turbulent models with small Mach num¬ 
bers the formation tends to be faster, because self-gravity dom¬ 
inates the large scales. The numbers of Class 0 protostars show 
only a narrow peak, followed by a similar, but shifted peak of 
Class 1 objects. In the models with higher turbulence, the ki¬ 
netic energy exceeds the gravitational one and the system is 
formally supported on global scales. Collapse only occurs lo¬ 
cally in the shock compressed cloud clumps. The formation of 
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Ser NGC7129 IC1396A Tau p Oph Cha I IC 348 


Fig. 2. The time r of the best correspondence between models 
and observations (in units of the global free-fall time rtf) for the 
seven investigated clusters. The symbols denote the parameters 
of the models: AJ = 0.5 (triangles), AJ = 2 (triangles down). 
Ad = 3.2 (squares). Ad = 6 (diamonds). Ad = 10 (stars); 1 < 
k <2 (open symbols), 3 < k < 4 (grey symbols), 1 < k < 8 
(filled symbols). The Gaussian collapse models are displayed 
as asterisks. 


Class 0 objects extends over a longer period (a few free-fall 
times) and in some cases there is a second burst of star forma¬ 
tion, following an increase of the number of prestellar cores, at 
a later time as in model M6k2a or M10k2. This second peak 
is not considered for the comparison with the observations, 
though. 

We count the numbers of objects in the particular classes 
and compare the relative numbers to the observational values 
from Table ^ We consider the entire observed population in¬ 
cluding more dispersed objects. Due to the use of periodic 
boundary conditions those are also included in the models. 

The time t of the best correspondence between observa¬ 
tions and models is determined when the weighted root mean 
square of the differences 


o-(f) = 



\n° - n™(f)Pw/ 


A 


2 

i=0 


( 1 ) 


becomes a minimum. The relative number of young stars in 
Class i (0, 1, 2) from observations is expressed as n° and n™ 
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denotes the relative number of YSOs in Class i from the mod¬ 
els. The factor w, is a weighting factor, introduced to account 
for the possible scatter due to small number statistics in both, 
the observations and the models. The weighting factor is set to 

OJi = 

The time t is shown in the lower panel of Fig. [2 and dis¬ 
played for all models in Fig. |2] In general. Serpens is fitted 
worse than the other clusters. The mean minimal deviation of 
all models is (cTnjjn) = 0.006 for p Oph, Cha I, and IC 348, 
(cTmin) = 0.008 for Taurus, (cr^jin) = 0.009 for IC 1396A, 
(cTmin) = 0.012 for NGC 7129, and (cTmin) = 0.017 for Serpens. 
In all cases cr(f) has a clearly defined minimum. The scatter of 
T is high. Its value varies from about 1.5 to 25 free-fall times 
for the different models, making it impossible to determine the 
actual ages of the clusters accurately. However, all models but 
one show the same sequence: T(Ser) < t(NGC 7129) < t(IC 
1396A) a; T(Tau) < rfp Oph) < T(Cha I) < t(IC 348). The 
only exception is model M6k2c, where t(IC 348) is slightly 
smaller than T(Cha I). In five models the age of IC 1396A is 
smaller than that of Taurus, in two models it is slightly larger, 
and in the other models the values are identical. Taking into 
account the second star-formation burst in the models M6k2a 
and M10k2 increases the ages of the older clusters, but it does 
not change the overall order. Thus, independent of the applied 
model. Serpens is the youngest cluster, NGC 7129 the second 
youngest cluster, Cha I and IC 348 are the most evolved clus¬ 
ters, while Taurus and IC 1396A are at roughly the same in¬ 
termediate evolutionary stage. However, the relative numbers 
of Class 0 and 1 objects in NGC 7129 and IC 1396A may be 
overestimated due to possible observational biases introduced 
by the Spitzer Space Telescope, as discussed in Sect. 12.61 These 
clusters may thus appear younger relative to the other ones. The 
fact that cr{t) always has one well-defined minimum and that 
all models produce the same evolutionary sequence indepen¬ 
dent of the initial conditions, supports a scenario of a single 
“burst” of uninterrupted, rapid star formation. This is in agree¬ 
ment with other theoretical and observational findings of star 
formation timescales and molecular cloud lifetimes (see e.g. 
Hartmann et al. 1200 II Ballesteros-Paredes & Hartmann l2005> . 
Several successive “bursts” of star formation within the cloud 
region are likely to alter the picture and make a determination 
of the age from the observed number ratios difficult. 

4.2. Star formation efficiency 

Figure |3 shows the star formation efficiency (the ratio of the 
mass accreted by the young stars to the total mass of the gas) at 
the time t versus the Mach number of the models for six of the 
seven clusters. We find an inverse correlation of the SFE with 
the Mach number (false alarm probability <0.5%). This is con¬ 
sistent with the theoretical findings by Klessen et al. J2000t and 
Heitsch et al. iEqoD, and probably due to the fact that in high- 
Mach number turbulence less mass is available for collapse at 
the sonic scale (Vazquez-Semadeni et al. l2003t . A linear fit is 
applied to the data and shown in Fig.|3 There is no correlation 
with the driving wave number. If we interpret the measured ve¬ 
locity dispersions in the clusters as the result of turbulence, we 


Table 2. The ratio of prestellar cores to the total number of 
YSOs at time r for three models and the observations 


Region 

M2k4 

M6k2a 

M10k2 

Observations 

p Oph 

0.30 

0.10 

0.10 

0.47 

Serpens 

0.79 

0.06 

0.74 

0.39 

Taurus 

0.32 

0.09 

0.22 

0.25 

Chal 

0.18 

0.11 

0.82 

< 0.39 


can estimate the SFEs at time r for the particular Mach num¬ 
bers from Eig. 0 In the case of p Oph this requires to extrapo¬ 
late the fit line beyond Al = 10. The SEE is ~ 0.27 in p Oph, 
~ 0.45 in NGC 7129, ~ 0.56 in IC 348, and between 0.60 and 
0.65 in Serpens, Taurus, and Cha I. (Lacking information on the 
velocity dispersion in IC 1396A, the corresponding SEE can¬ 
not be calculated.) These values are significantly higher than 
the measured SEEs, which are only around 0.1. (No informa¬ 
tion is available on the SEEs of IC 348 and NGC 7129, but we 
expect them to be in the range of the other clusters.) Only in the 
case of p Oph the SEE of the models is in the range of the SEE 
measured. The main reasons for this discrepancy are probably 
the limitation of the gas reservoir and the neglect of outflows 
and feedback mechanisms as well as of magnetic fields in the 
simulations. Bipolar outflows limit the local SEE, because the 
protostellar jet will carry a certain fraction of the infalling ma¬ 
terial away, furthermore, its energy and momentum input will 
affect the the protostellar envelope and may partially prevent it 
from accreting onto the protostar (e.g. Adams & Eatuzzo fl996> . 
The presence of magnetic fields would also retard the conver¬ 
sion of gas into stars (for current simulations see Heitsch et al. 
120011 Vazquez-Semadeni et al. 120051 Li & NakamuraEQQl. 

4.3. Presteiiar cores 

As an additional test, the numbers of prestellar cores are anal¬ 
ysed as described in Sect. |3for the three models shown in Eig.^ 
and compared to the four star-forming regions, where informa¬ 
tion on the number of prestellar cores is available. The ratios 
of the number of prestellar cores to the total number of YSOs 
(Upsc/n*) at time t are listed in Table|3together with those val¬ 
ues from the observations (calculated from Tabled). The ob¬ 
served ratio in Taurus are roughly represented by models M2k4 
and M10k2, but for the other regions the models produce ei¬ 
ther a much higher or much lower ratio. We intended to check 
if Mpsc/^psc allows us to draw conclusions about the number 
of prestellar cores that actually collapse and form stars. This 
seems not to be possible on the basis of the current analysis. 
In addition, defining a prestellar core is rather difficult, both, 
from an observational and a theoretical point of view. This adds 
another level of uncertainty to results based on prestellar core 
statistics. 

5. Summary and conclusions 

We analysed the temporal evolution of the fractions of YSO 
classes in different gravoturbulent models of star formation and 
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Fig. 3. The star formation efficiency at time t plotted versus the Mach number of the models. The symbols are the same as in 
Fig.El circles denote models with Al = 0.1. The dotted lines show a linear fit to the data. 


compared it to observations of star-forming clusters. The ob¬ 
served ratios of Class 0, 1, 2/3 objects in p Ophiuchi, Serpens, 
Taurus, Chamaeleon I, NGC 7129, IC 1396A, and IC 348 
can be reproduced by the simulations, although the time when 
the observations are best represented varies depending on the 
model. Nevertheless, amongst the clusters with good observa¬ 
tional sampling we always find the following evolutionary se¬ 
quence of increasing age: Serpens, Taurus, p Oph, Cha I, and 
IC 348. 

We find an inverse correlation of the star formation effi¬ 
ciency with the Mach number. However, our models fail to re¬ 
produce the observed SFEs for most of the clusters. This is 
probably due to the lack of a sufficiently large gas reservoir in 
the simulations and the neglect of energy and momentum input 
from bipolar outflows and/or radiation from young stars. Only 
the SFE in p Oph is reproduced. This region is characterised 
by very high turbulent Mach numbers. The fact that our simple 
gravoturbulent models without feedback are able to reproduce 
the right number ratios of YSOs for Ad ^ 10 suggests that pro- 
tostellar feedback processes may not be important in shaping 
the density and velocity structure in star-forming regions with 
very strong turbulence and argues for driving mechanisms ex¬ 
ternal to the cloud itself (see also the discussion in Ossenkopf 
& Mac Low l2002l and Mac Low & Klessen l2004t . 

The relative numbers of YSOs can reveal the evolutionary 
status of a star-forming cluster only with respect to other clus¬ 
ters, the absolute age is difficult to estimate. Better agreement 
between models and observations requires a better considera¬ 
tion of environmental conditions like protostellar outflows and 
magnetic fields in the simulations. A larger observational sam¬ 
ple, achieved by complete censuses of more star-forming re¬ 
gions, is also desirable. 
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